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ABSTRACT 

We derive mass functions (MF) for halo red dwarfs (the faintest hydrogen burning stars) and 
then extrapolate to place limits on the total mass of halo brown dwarfs (stars not quite massive 
enough to burn hydrogen). The mass functions are obtained from the luminosity function of a 
sample of 114 local halo stars in the USNO parallax survey (Dahn et al. 1995). We use stellar 
models of Alexander et al. (1996) and make varying assumptions about metallicity and about 
possible unresolved binaries in the sample. We find that the MF for halo red dwarfs cannot rise 
more quickly than 1/m 2 as one approaches the hydrogen burning limit. Using recent results 
from star formation theory, we extrapolate the MF into the brown-dwarf regime. We see that 
likely extrapolations imply that the total mass of brown dwarfs in the halo is less than ~ 3% of 
the local mass density of the halo (~ 0.3% for the more realistic models we consider). Our limits 
apply to brown dwarfs in the halo that come from the same stellar population as the red dwarfs. 



dark matter - Galaxy: stellar content - Galaxy: halo - stars: low mass, brown dwarfs - stars: 
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1. Introduction 

In this paper we determine a mass function for low mass stars and then extrapolate to derive an upper 
limit on the number of brown dwarfs. A key quantity for describing a stellar population is the mass function 
£(to), the density of stars with mass between m and m + dm; for £(m) we use units #stars/pc 3 /M0. The 
mass function (hereafter, MF) of halo stars can address the question of what makes up the dark matter 
of our Galaxy. Although it has been clear for some time that main sequence stars heavier than our sun 
are not sufficiently abundant to explain the mass of the Galaxy, the question has remained whether or not 
there can be large numbers of lower mass objects. Very low mass stars or substellar objects have been 
considered the most plausible candidates for baryonic dark matter. [White dwarfs, another possibility for 
baryonic dark matter, are not considered in this paper.] Recent work has ruled out red dwarfs, stars just 
barely massive to burn hydrogen (m > O.O92M0), as a significant source of dark matter (Graff & Freese 
1996; Boeshaar, Tyson & Bernstein 1994; Bahcall et al. 1994). Brown dwarfs, star-like objects that are just 
barely too small to burn hydrogen, are the remaining low mass candidates. These are extremely difficult to 
observe optically (see e.g., Burrows and Liebert 1993 for a review) although there have been searches for 
them with gravitational microlensing (Alcock et al. 1995, Ansari et al. 1996). 

As our main result, we determine a MF of halo red dwarfs, and then extrapolate into the brown dwarf 
regime to obtain an upper limit on the number of halo brown dwarfs. If brown dwarfs are a primary 
component of halo dark matter, then we would expect a steeply rising MF as we go to decreasing mass. 
Instead, we show that likely extrapolations imply that brown dwarfs are not nearly numerous enough to be 
a significant component of halo dark matter. As a caveat, let us point out that our limits apply only to 
brown dwarfs in the halo that come from the same stellar population as the red dwarfs (see the discussion 
at the end of the paper) . We also rely on recent results from star formation. The goal of our paper is to 
show that brown dwarfs coming from known stellar populations are not present in great abundance. ^ 

Unfortunately, the MF cannot easily be obtained by observations, since it is difficult to measure the 
mass of a star. Instead, we must convert from the measured quantity, the luminosity function ip(M), the 
number of stars in a magnitude range M — > M + dM (note that M refers to magnitude while m refers to 

1 There have been several attempts to determine the mass function of local stars (Miller & Scalo 1979, 
Kroupa, Tout & Gilmore 1993, Henry 1995). Local stars are predominantly Population I disk stars. A 
population II MF for low mass stars from globular cluster NGC 1261 can be found in Zoccali et al (1995). 
In our paper we focus instead on local field halo stars, which are not in globular clusters. 
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mass). For the luminosity function (hereafter LF), we use units #stars/pc 3 /M. 

To obtain a mass function for red dwarfs, we use the LF from the extensive parallax survey being carried 
out by the US Naval Observatory. Preliminary results were reported in Monet et al (1992) with a more 
complete update in Dahn et al (1995). They published an LF based on 114 stars with velocities tangential 
to the line of sight Vtan > 220km/sec. These high velocity stars should all be halo stars (temporarily in 
the disk) since true disk stars have a local velocity dispersion of only ^30 km/sec. Mera, Chabrier & 
Schacffer (1996) used this LF to derive an MF as discussed further below; see their Fig. 1. In this paper 
we present a more complete derivation of the mass function for red dwarfs, and a different interpretation. 
We take into account several potentially important effects not considered by Mera et al, such as unresolved 
binaries in the sample and the effects of different choices of metallicity. In particular, many of the 'stars' in 
the sample are likely to be unresolved binaries, and substantial numbers of binaries may drastically alter 
the results. The MF derived by Mera et al falls into the range of possibilities bracketed by the parameter 
range we consider. However, our interpretation is different. We emphasize a different slope of the MF near 
the hydrogen burning limit. (By "slope", we mean the logarithmic slope dlog£/dlogm.) In addition we 
use theoretical results from star formation to extrapolate to lower masses and limit the abundance of halo 
brown dwarfs. Because we are using parallax data, we can only discuss the local population of halo stars. 



2. Deriving the MF of Red Dwarfs 

As described above, we use the LF obtained from the parallax survey of Dahn et al (1995). To convert 
from the luminosity function vp(My) to the mass function £(rn), we use 

V(M y ). (1) 

Clearly we need a mass-luminosity relation. We use the models calculated by Alexander et al (1996). In 
Fig. 8 of their paper, one can see that the mass-luminosity relation is non-linear. The luminosity of stars 
drops fairly slowly with decreasing mass until an inflection point at m « 0.4M Q . Then the mass- luminosity 
relation turns over, and the luminosity rapidly becomes dimmer as the mass decreases toward win ■ We fit 
between the calculated points of (Alexander et al. 1996) using cubic splines. Since our MF will depend 
on certain assumptions, we will always consider the most 'conservative' case, i.e., we generate the steepest 
possible MF as one goes to low mass, corresponding to the largest number of low mass red dwarfs and of 
brown dwarfs. 



£(m) 



dM v 



dm 
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We use the model of Alexander et at with stellar metallicity Z = 3 x 10 -4 . Nearly all stars in the 
sample have larger metallicities than this, as can be seen from the comparison of parallax data points with 
theory for various metallicities in Fig. 6 of Alexander et al. (1996). We use this value of Z because it gives 
the most conservative result (the steepest MF towards low mass). One can see this effect in Fig. la for two 
different metallicities and can understand it as follows. High Z stars are dimmer than low Z stars of the 
same mass, while high Z stars are more massive than low Z stars of the same luminosity. For the range of 
luminosities we are considering, if we underestimate Z, stars would be on a much steeper part of the My{m) 
curve and many would have masses in a narrow mass range close to the hydrogen burning limit TOjj. In 
addition, the factor \dMy /dm\ of Eq. (|j) becomes large for stars near run (because the mass range for low 
Z stars of a given luminosity is so small). These two effects combine to cause the MF to steepen towards 
lower mass as we lower the metallicity. Thus, using Z = 3 x 10~ 4 gives the steepest allowed MF and the 
largest number of brown dwarfs. 

In addition, the reported LF can differ from the true LF because of potential errors and biases in the 
parallax data. Observers on average underestimate the luminosity of the stars because there are more 



distant stars than nearby stars ( Lutz fc Kclkcr 1973 ). Also, due to random errors, there is a diffusion-like 



broadening of the LF ( Btobic, Ishida fc Peacock 1989f ) so that the observed LF does not drop as fast as the 



true one as one goes to dim stars. Both systematic biases (known as Malmquist bias) will cause us to have 
two many small dim stars in our sample, so we will over-estimate the steepness of the MF as one goes to 
small mass. Again, our results are conservative. 

The LF derived by Dahn et al. grows to a maximum value at My = 12 and then drops off again as one 
goes to lower mass. [Hereafter, we use the words rising and falling to refer to the behavior of functions as 
one goes to smaller values of mass or of luminosity] One might naively interpret this to mean that there 
cannot be many stars of low mass. However, this interpretation may not be correct because of the shape of 



the mass-luminosity curve (Kroupa, Tout & Gilmore 199C); 5LF that reached a maximum and then turned 



over for all a decreasing luminosity function does not necessarily imply a decreasing MF. 

A Complication: Binaries: A complication arises in the implementation of Eq. (jl]): some of the 
stars in the survey may actually be unresolved binaries. If so, the luminosity of the binary is due to light 
from both stars. If we mistakenly interpret the light to be only from a single star, we may overestimate the 
mass of the star. To deal with this potential problem, we make a variety of assumptions about unresolved 
binaries, and plot the corresponding mass functions in Fig. 1. Our models represent extremes which bracket 
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the possible range of MFs due to unresolved binaries. We shall see that for all cases, we still conclude that 
brown dwarfs are an insignificant component of the mass of the halo. 

In Model I we make the most common assumption about binary stars: the masses of the two stars are 
independent and yet are drawn from the same distribution; i.e., there is a single MF which determines the 
mass of each star. Note that we do allow the binary companion to have a mass below the hydrogen burning 
limit (we believe that there is no physical reason to cut off the distribution here). Thus our model (and 
its conclusions) is different from that of other authors (Kroupa, Tout & Gilmore 1991, 1993; Piskunov and 
Malkov 1991) who do not consider brown dwarf companions. We here show that if the mass density in 
brown dwarfs is very small for the no binary case, then it is also very small for the case of model I. Assume 
that brown dwarfs are far more numerous than true stars. Then by our assumption that the masses of the 
two stars are drawn from the same MF, the secondary in any binary system is nearly always |^] a brown 
dwarf |J The contribution of the brown dwarf to the total luminosity of the system is negligible and we can 
derive the mass of the primary from the luminosity of the system. This means that the MF derived using a 
single star mass-luminosity relation also describes the MF for binary stars. But, as we will show below, this 
single star MF implies that there are not large numbers of brown dwarfs, contradicting our assumption. 
Thus, if the mass density in brown dwarfs is very small for the no binary case, then it is also very small for 
the case of model I. 

Although model I is the most plausible, it is difficult to directly calculate the actual mass function 
from the observed massfunction (although the reverse process is trivial). Thus we tried two other binary 
star models. Model II is plausible. Model III is extremely unlikely, and is designed to exaggerate the shift 
in slope of the MF. In both models, we analytically determined the MF which would produce the observed 
LF and which would satisfy the assumed binary composition. Both models steepen the MF towards low 
masses. [Note that the direction of this effect is in agreement with what was found by Kroupa, Tout, & 

2 Since the above authors did not allow brown dwarf companions, this statement would not be true for 
their work; our model II more closely approximates their results. 

3 There has been an intensive search for brown dwarfs as binary companions. This search has only found 
two brown dwarfs despite examining hundreds of systems; thus brown dwarfs make up < 1% of the mass of 
disk stars (Rebolo, Zapatero Osorio & Martin 1995; Nakajima et al. 1995). Because brown dwarfs cool as 
they age, one can only hope to optically detect young brown dwarfs. Thus, these optical searches are not 
sensitive to a halo population of brown dwarfs. 
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Gilmore (1991, 1993), who studied the effects of unresolved binaries in other data sets (typically Pop. I)]. 

In model II, each "star" has a probability / (the binary fraction) of being an unresolved binary. In 
every binary system, we take the mass of the secondary m s to be evenly distributed between zero and the 
mass of the primary to p ; i.e., the binary mass ratio q = ^ has an equal chance of being anywhere in the 
interval [0,1]. For each bin My in the observed LF histogram tp/(Mv), a fraction (1 — f)ip'(Mv) of the 
LF is assumed to represent single stars whose absolute magnitude has been accurately determined. The 
remaining fraction fij)'{Mv) represents all unresolved binaries and is divided into 100 parts i, each having 
a mass ratio qi = i/100, where i is an integer between and 100. For each part i we find the two masses 
TOp.i and m St i = qim p ^ with visual magnitudes My P ,i and My s ,i such that the combined luminosity of the 
system is exactly what is observed. We then put these pieces together to obtain the luminosity function 
ijj(Mv) corrected for unresolved binaries; the results from model II will be discussed shortly. 

In model III, we assume that all of the stars are in binary systems, and that the two stars have the 
same luminosity. Thus the luminosity of the primary is half that of the system. Compared to the no-binary 
case, in model III the luminosity of each star is half as much, the magnitude is larger by an amount 0.75, 
and the mass of each star is smaller. Also, there are twice as many stars because each system is a binary. 
Thus, according to model III, the "true" LF ip(Mv) can be derived from the observed LF ip'(Mv) by: 
ip(Mv) = 2 x ip'(Mv — 0.75). Although the shape of the LF is unchanged, this model generates a steeper 
MF for the same reasons a lower metallicity model generates a steeper MF: stars lie on a steeper part of the 
My (to) curve so that their masses lie in the narrow mass range near the hydrogen burning limit, and the 
factor \dM v /dm\ of Eq. (§) becomes large as discussed in section 2. 

Our results for all the models are plotted in Fig. 1. The MF of red dwarfs has been obtained from the 
LF using Eq. ([j]) for four different cases. Both binary models II and III generate MFs that are steeper than 
if we ignored unresolved binaries. In the plots, the MF has been multiplied by to 2 so that one can most 
easily see how it approaches the hydrogen burning limit (as described in the next section). Our result is 
that, at the low mass end, the MF behaves as £(m) ~ mT a with a < 2 for all four cases we consider. In 
fact, only binary model III could have a — 2 while all other cases have a < 2. 

For comparison, Mera et al (1996) found the MF only for the case of no binaries, and roughly estimated 
a ~ 2. They obtained this number by fitting to the MF over a large mass range. However, we find that for 
this case of no binaries, the slope a is smaller near the hydrogen burning limit, as one can see from Fig. 
la; it is this smaller slope at the lowest mass end of the red dwarf range that is relevant for extrapolating 
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into the brown dwarf regime]]. Below we also take advantage of recent results from star formation theory 
to guide our extrapolation, and this leads us to a different interpretation of the importance of the shape of 
the MF near the hydrogen burning limit. Also, our work has considered a wider range of possibilities for 
binaries with the aim of obtaining a careful estimate of the range of possibilities for a. As seen in the next 
section, the fact that we find a slope which always satisfies a < 2 has important consequences. 



3. Upper Limit on the Total Mass of Halo Brown Dwarfs 

Extrapolating the MF to Brown Dwarfs Theoretical Input: Given our mass functions of 
halo red dwarfs, we will now extrapolate to the lower mass brown dwarfs. This extrapolation should be 
guided by theory. The work of Adams and Fatuzzo (1996) predicts the mass of a star based on the physical 
properties of the cloud core that forms that star. In their theory, stars determine their own masses through 
strong stellar winds and outflows. They predicted two extremes for the MF of stars forming in molecular 
clouds. At one extreme, the stellar mass is determined by one factor, e.g., the sound speed of the molecular 
cloud core. In this extreme, the mass function of stars would be a power law. At the other extreme, the 
mass of the star is determined by a large number of independent physical variables. Invoking the central 
limit theorem, Adams and Fatuzzo reasoned that log£(?7j) should be a Gaussian; i.e., the mass should be 
log-normal distributee^]. 

These two cases represent extreme limits. We expect the actual function to lie somewhere between the 
two. In short, a MF should not be concave up, only flat or concave down on a log-log plot. There are other 
current theories of star formation based on completely different physical reasoning, such as heirarchical 
fragmentation ( |Zinneckcr 1984 ) and collisions between cloud cores ( Price fc Podsiadlowski 1995 ). All 



of these theories predict mass functions that rise less quickly than power laws, and thus lie within the 
same range of possibilities as discussed above. All the possible MFs generated by our different models are 
consistent with a mixture of the two theoretical models over the red dwarf mass range, i.e., all are flat or 
concave down. 



4 Mera et al (1996) found a ~ 2.5 for the data of Richer and Fahlman (1992); however these data had 
difficulty distinguishing least luminous stars from galaxies and disagree with the results of Hubble Space 
Telescope (Bahcall etal 1994, Graff and Freese 1996). 

5 To leading order, the MF of disk stars appears to be fit by a log-normal (Miller & Scalo 1979). 
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Upper limit on the total mass in Halo brown dwarfs: The total mass in brown dwarfs is 

rntot = / m£(m)dm. (2) 
Jo 

Since there are so few red dwarfs relative to the halo, the total mass of brown dwarfs can only be large if the 
MF is steeply rising as one goes to small masses. In the theory discussed above, any possible extrapolation 
of an MF lies between two possible extremes: power law or log-normal. Within these two bounds, there will 
be more brown dwarfs if we extrapolate with a power law. Hence we use £ ex m~ Q to place an upper bound 
on the total mass of brown dwarfs in the halo. Note that we extrapolate the slope at the low mass end of 
the observed mass function. 

If a < 2 then the integral in equation ^ converges and one finds the total mass in brown dwarfs to be 
of the same order of magnitude as the total mass in red dwarfs. If a > 2, then the integral diverges. Since 
the mass density of brown dwarfs is not infinite, we mean by this divergence that the total mass density is 
dependent on a lower mass cutoff of the MF. If a = 2 then the integral diverges, but only logarithmically, 
i.e., each order of magnitude of mass range contains an equal total mass. Thus, even for a lower limit of 
10 _7 M Q , the total mass in brown dwarfs ~ 12 times the mass in red dwarfs, or only 3% of the local halo 
mass. If a > 2, the integral diverges fast enough that there could be enough brown dwarfs to fill the halo. 

As shown in Fig. 1, for all of our models, with varying metallicities and binary content, a power law 
extrapolation to low masses has a power of a < 2. Thus brown dwarfs cannot contribute more than ~ 3% 
of the local halo mass. Only our most extreme model (binary model III) can have a = 2; the more sensible 
models have a < 2. Eq. (0) with a < 2 implies that there cannot be many more brown dwarfs than red 
dwarfs. Previously we estimated that red dwarfs contribute ~ 0.3% of the halo (for a halo mass density of 
7 x 10 -3 Mq/pc 3 ). In summary, brown dwarfs cannot contribute more than ~ 3% of the mass of the halo, 
and in realistic models the limit is an order of magnitude smaller. 



4. Discussion 



Red dwarfs are a tiny fraction, 0.3%, of the mass of the galactic halo ( Graff fc Freese 1996 ). In 
order for any appreciable fraction of the halo mass to consist of low mass star-like objects, there must be 
substantially more mass in brown dwarfs than in red dwarfs. We estimated the mass in brown dwarfs by 
extrapolating the mass function (MF) of red dwarfs below the hydrogen burning limit. Our extrapolation 
was based on current theories of star formation, which predict that the MF is either flat or concave down 
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on a log-log plot (or in fig. 1), i.e., the MF rises no faster than a power-law. We made several assumptions 
about the metallicity and binary composition of the stars in the sample. We found that likely extrapolations 
imply that brown dwarfs make up less than ~ 3% of the local mass density of the halo. For our most 
realistic models, the limit on the total mass in brown dwarfs is roughly the total mass in red dwarfs, ~ 0.3% 
of the local mass density of the halo. 

To repeat an earlier caveat: our limits on the brown dwarf density assume that the brown dwarfs in 
the halo come from the same stellar population as the red dwarfs. The possibility always remains that there 
are large numbers of brown dwarfs from an entirely different population of stars (Population III). However, 
theoretical work on star formation (Adams and Fatuzzo 1996) indicates that an earlier population of stars 
is likely to be skewed towards a predominance of high mass stars, not brown dwarfs. There may also be 
some yet undiscovered mechanism of star formation which allows the MF to rise faster than power-law. 
In this paper we have presented red dwarf mass functions and have shown that brown dwarfs which are 
members of known stellar populations are not present in the halo in great abundance. 

We thank Santi Cassisi for sharing unpublished work with us. We also thank F. Adams, G. Basri, P. 
Kroupa and J. Bahcall for useful discussions. We acknowledge support from NSF PHY-9406745 and from 
the Univ. of Michigan Physics Dept. 
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FIGURE CAPTIONS 

Figure 1: The mass function of red dwarf halo stars. Each of the four models is derived from the 
LF of Dahn et al (1995) but assumes different metallicity and binary content. In all three panels, crosses 
without errorbars illustrate the mass function derived for stars with metallicity Z = 3 x 10~ 3 and no binary 
companions. Error bars are due to Poisson errors in the LF. The other model presented in panel (a) has 
Z = 6 x 10~ 3 (no binaries) for comparison. Lower metallicity leads to a MF with more low mass stars. 
Panels (b) and (c) show binary models II and III for Z = 3 x 10 -3 as described in the text and can be 
compared with the no binary model (crosses) of the same metallicity. For binary models we plot the total 
MF (primaries and secondaries). Panel (b) shows model II for two different binary fractions, / = {0.5, 1.0}; 
/ = 0.5 is closer to the no-binary model (crosses). Binary model II has a somewhat steeper MF than the 
no binary case. Binary model III has been designed to exaggerate the number of low mass stars compared 
to high mass ones and is unrealistic. 

We have multiplied the vertical axis by m 2 to emphasize that all MFs converge at low mass (see eq 
and subsequent text). A MF which is decreasing to the left will converge; one that is increasing will 
diverge. One that is flat will diverge, but only logarithmically. Even for the extreme case of binary model 
III, the MF is flat or decreasing and the total mass in brown dwarfs could be at most ~ 3%. For more likely 
models I and II, the limit on the total mass in brown dwarfs is very roughly the total mass in red dwarfs, 
~ 0.3% of the local mass density of the halo. 



-4 







GO 
GO 

£ 

X 



-5 - 



-4 



O 



CO 

I 

O 
GO 

■+-> 

GO 

=tfc: 



o 



-5 



-4 



-5 - 



III 
a) 

> 

X 


1 1 

me 

X J xj T t 
X TXXXJ X I xl X } 

fftr -- 


1 1 J_ 
tallicity 

■sr 

- x__ x 


~b) 

1 

j 

> 
X 


mod* 

1 1* i 1 1 I f I - 


3l II 

z x 
k 

1 


~c) 

II 

> 
X 

1 1 1 


lit * I T mc 

x x x X x I x 

1 


.del III 

x x 

II 

1 1 1 



0.1 



mass (M Q ) 



